Topology of Fracture Networks 
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We propose a mapping from fracture systems consisting of intersecting fracture slieets in three 
dimensions to an abstract network consisting of nodes and links. This makes it possible to an- 
alyze fracture systems with the methods developed within modern network theory. We test the 
mapping for two-dimensional geological fracture outcrops and find that the equivalent networks are 
small-world and dissasortative. By anlayzing the Discrete Fracture Network model, which is used 
to generate artifical fracture outcrop networks, we also find small world networks. However, the 
networks turn out to be assortative. 

PACS numbers: 89.75.Hc,91.55.Jk,81.40.Np,89.75.Kd 



Topological analysis of networks has had an explo- 
sive growth over the last decade [l|. A large number 
of new concepts and quantitive tools for describing net- 
works have been introduced, making it possible to de- 
scribe and classify complex network structures at a level 
that never earlier has been achieved |^, y] • There is one 
class, though, of networks that has resisted this kind of 
analysis: Fracture networks. These consist of intersect- 
ing fracture sheets, making both the concepts of links and 
nodes far from obvious. Fracture networks, however, are 
extremely important from a technological point of view. 
For example, in carbonate petroleum reservoirs, the oil is 
transported through fracture networks as the permeabil- 
ity of the porous matrix is too low [J] . Another example 
is the extraction of shale gas though hydrofracturing [5| . 

We propose a transformation from fracture network 
to an equivalent network consisting of nodes and links. 
This makes it possible to qualitatively and quantitatively 
characterize the topology of fracture networks. 

An important consequence of this is that it is possible 
to compare models that generate artificial networks with 
real networks quantitatively. 

Fracture outcrop networks have been studied from a 
network point of view by Valentini et al. [a, lZ|. Frac- 
ture outcrops are fracture lines visible on the surfaces 
of geological formations. The outcrop fracture lines are 
one-dimensional cuts through the two-dimensional frac- 
ture sheets. Valentini et al. treats fracture lines as links 
and their crossing points as nodes, this gives a more nar- 
row degree distribution than the transform proposed in 
this paper. However, Valentini et al. also conclude that 
fracture networks are small- world networks |6J- In three 
dimensions where the fractures are sheets, the transfor- 
mation we propose is necessary to define the topology 



network. 

Our analysis is somewhat related to the information 
measure for cities introduced by Rosvall et al. [8|. 

Lacking data on three-dimensional fracture systems, 
we analyse in the following fracture data from eight out- 
crops found in south-east Sweden. A detailed descrip- 
tion of the bedrock composition and geological history 
are given in [9|, llO| . We show one of the outcrop frac- 
ture networks in Fig. [1^. As we shall see, the equivalent 
network (shown in Fig. [TJd) constructed from the origi- 
nal network has small- world character. Furthermore, it 
is disassortative. 

We then go on to analyse artificial fracture networks 
generated with the Discrete Fracture Network (DFN) 
model 
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The equivalent networks constructed from 
the original networks generated by this model also show 
small-world behavior. However, they are assortative. 

The eight outcrops covers between 250 and 600 m^. 
All visible fractures with length over 0.5 m have been 
recorded in the data sets. We prepare the data sets 
as follows. When tracing the fracture lines, they may 
appear disconnected or doubled due to topography or 
ground weathering. An illustration of a outcrop is shown 
in Fig. [2^1. We therefore use a reconnection procedure 



12l |. That is, we first project fracture traces on a flat 



surface to reduce the perturbation due to rock surface 
topography. Then scattered segments that are likely to 
belong to the same trace are reconnected to one single 
segment accounting for orientation and distance consis- 
tency. We focus on traces with a dashed-line, discon- 
nected step or layered patterns. We then straighten all 
the fractures lines. The result is shown in Fig. [2}d. 

We have now come to the central idea of this paper. 
In Fig. [2j;, each fracture line has been associated with a 









FIG. 1: (Color online) a) Fracture network of outcrop 
AMS000025 b) Equivalent network based on the original net- 
work shown in a). 



node. Whenever two fracture lines cross, we place a link 
between the nodes representing the two fracture lines. In 
Fig. [2Ji, we show the equivalent network consisting of 
nodes representing the fracture lines and links represent- 
ing crossing fracture lines 13 1. 

We note that this equivalent network is as simple 
to construct in a three-dimensional system of fracture 
sheets: each fracture sheet is represented by a node and 
whenever two sheets cross each other, a link is placed 
between the equivalent nodes. 

Arguably the most central property of any complex 
network is the degree distribution p{k). The degree, fc, 
of a node is the number of other nodes that it is linked 
to. The equivalent networks generated from the outcrop 
networks show a broad degree distribution. Wc plot the 
cumulative distribution, P{k) in Fig. [31 When P{k) fol- 
lows a power-law the network is scale free [2|. We plot 
the data against log-log fits, but note that our geological 
data can be well fitcd by a log- normal. In the case of the 
DFN data the goodness of the log-log fit is dependent on 
model parameters. 

The clustering is a local measure of how well a network 



FIG. 2: (Color online) Clock-wise from upper left, a) Repre- 
sentation of fracture outcrop network, b) Reconnected frac- 
ture network, c) Equivalent network placed on top of fracture 
outcrop network, d) Equivalent network representation of b). 
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FIG. 3: (Color online) Cumulative degree distribution for: a) 
Networks generated from the eight outcrop data sets. Insert 
shows the average compared to a log-log (exponent —2.3) and 
a log-normal fit. b) DFN model. Values for the expoents of 
the fits are given in Table HIl 



is connected on a local neighbor-to-neighbor scale. The 
global clustering coefficient, C, is defined [ij, |l5| as the 
average over all the local clustering coefficients, Q, for 
each node 
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where ki is the degree of node i^ N is the total number 
of nodes and Epf^^t is the number of links between the 
nearest neighbors of node i. The clustering coefficient 
falls in the interval < C < 1 , and a high value indicates 
that there is a high chance that two neighbors of a node is 
connected to each other. This makes the network highly 
connected on a local scale, making it easy for nodes to 
efficiently interact on this scale. 

In order to determine whether the clustering coeffi- 
cients found for the networks are large for their number of 
nodes and links, we compare them to rewired and random 
versions of the same networks. In rewiring [16[ two pairs 
of connected nodes are selected at random, and the links 
interchanged so that two new pairs of connected nodes 
are created. The procedure is repeated until all links 
are moved. This preserves the degree distribution since 
all nodes retain their initial degree, but it removes any 
correlation between the degrees of the connected nodes. 
For the random version all links are removed and redis- 
tributed randomly between the nodes. This produces a 
new degree distribution that is generally not broad. In 
all cases the quoted values for these networks are aver- 
aged over 1000 realizations. As can be seen from TableU 
the equivalent networks have an average clustering coef- 
ficient of 0.18 which is more than an order of magnitude 
larger than for comparable rewired networks, and two or- 
ders of magnitude larger than for purely random versions. 
Hence, they are well connected on a local scale. 

The efficiency, E, is a global measure for how well the 
different parts of the network are connected, and how 
easily nodes in different parts of the network can interact. 
The measure is defined usiiig the shortest distance, d^, 
between two nodes i and j [3| 
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where dij = oo if node i and j are not connected. E falls 
in the interval < E < 1 , and a high value indicates that 
it is easy for nodes far apart in the network to interact 
since there on average is just a few links between any two 
nodes. 

In TableU we present E for all the equivalent networks 
and their average is 0.065, which is smaller than for the 
rewired (Erw) and random {Era) versions both having 
an average of 0.11. However the efficiency {E) is only 
smaller by a factor of about 2, making E and Er^/ra 
of the same order. We would expect the rewired and 
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FIG. 4: (Color online) Plot of the correlation matrix C(fci, k2) 
based on the equivalent networks generated from the eigth 
outcrop fracture data sets. 



random networks to have a high efficiency, several orders 
of magnitude larger than ordered networks, because they 
have a large portion of long-range links. The fact that the 
equivalent networks have an efficiency comparable to that 
of the rewired and random versions means that compared 
to ordered networks they have a large efficiency. We will 
discuss the impact of C and E for the equivalent networks 
in more detail below. 

It is also interesting to study any correlations between 
the degrees of linked nodes. Does high degree nodes link 
predominantly to low degree nodes or high degree nodes? 
Maslov and Sncppen [17| introduced a correlation matrix 
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where P{ki, k2) is the probability that a node of degree 
ki is linked to a node of degree k2 for the network to 
be investigated. PR{ki,k2) is the same probability of a 
rewired version of the network. If C(fci,fc2) ~ 1 for all 
(^1,^2) then there is no degree correlations in the link- 
ing between nodes. If C(fci,fc2) > 1 for some values of 
(^1,^2) then there is an over-representation of links be- 
tween nodes of degree ki and fc2 in the investigated net- 
work compared to that of a rewired version of the net- 
work. If C(fci,fc2) < 1 there is an under- representation. 
Note that the matrix C{ki,k2) is symmetric. 

In Fig. |4] we have plotted the average of the ma- 
trix C(fci,fc2) for all outcrops, where -Pfl(fci,fc2) is av- 
eraged over 10000 realizations. We observe an over- 
representation of small degree nodes linking to higher de- 
gree nodes, and an under-representation of equal degree 
nodes linking to each other. Such networks are disassor- 
tatiye, and are abundant in naturally occurring networks 

The characteristic path length, L is defined as the av- 
erage distance between any pair of nodes of a network, 
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TABLE I: List of the number of nodes (fractures), links, maximum degree kmax, average degree k, clustering coefficient C, 
clustering coefficient for rewired networks Crw, clustering coefficient for random networks Cra, efficiency E, efficiency for 
rewired networks Erw, and efficiency for random networks Era for all the outcrop samples. 



Sample 



Nodes Links k„ 



C Crw Cra E Erw Era 



AMS000025 


787 


858 


23 


AMS000026 


716 


520 


20 


AMS000205 


973 


1188 


32 


AMS000206 


737 


487 


11 


AMS000208 


955 


1297 


31 


AMS000209 


955 


1162 


27 


AMS100234 


946 


1549 


44 


AMS100235 


785 


1392 


44 



2.18 0.170 0.0048 0.00178 0.046 0.104 0.101 

1.45 0.088 0.0033 0.00087 0.019 0.048 0.032 

2.44 0.193 0.0043 0.00174 0.032 0.122 0.118 

1.32 0.120 0.0013 0.00067 0.004 0.033 0.020 

2.72 0.226 0.0067 0.00213 0.079 0.138 0.138 

2.43 0.177 0.0050 0.00178 0.068 0.119 0.118 

3.27 0.236 0.0138 0.00291 0.133 0.164 0.172 

3.55 0.243 0.0180 0.00394 0.141 0.176 0.192 



Average 



857 1057 29 2.42 0.182 0.0072 0.00198 0.065 0.113 0.111 
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FIG. 5: (Color online) Examples of fracture systems gener- 
ated with the DFN model for varying parameters ai and D2 ■ 



Having a large clustering coefficient indicates a large lo- 
cal connectivity, and a small characteristic path length 
indicates a large global connectivity. When both of these 
criteria are fulfilled, wc have a small- world network 14 1 . 



Networks consisting of more than one disjoint part will 
have dij = 00 for at least one pair of nodes. Hence, the 
characteristic path length is not a good measure for the 
global connectivity of such networks. However a small 
value of dij for most pairs of nodes will give a large av- 
erage value for l/d.y which is measured by the efficiency. 
Therefore a large E is comparable to a small L for de- 
scribing the global connectedness. Since the fracture net- 
works found in the outcrops have been shown to have a 
clustering coefficient significantly larger than rewired and 
random versions, and an efficiency of the same order as 
the rewired and random networks wc conclude that these 
are small- world networks. 

We now turn to analyzing the DFN model [ll[. It is 



based on the observation that the length of fracture lines 
in outcrops, Z, arc distributed according to a power law 
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a/ = 3 (ASM000205 and ASM000206) and one with 
ai = 2.3 ± 0.2 (the rest) [2^. The angular distribution 
of the directions of the fractures depends on the fracture 
system. We assume here the simplest, i.e., a uniform dis- 
tribution. The outcrop data studied show strong signs of 
prefered directions for the fractures, but using this in the 
angular distribution of the DFN model does not have a 
significant impact. Lastly, the position of the fractures 
must be specified. The DFN model uses a hierarchical 
construction 23|, |2J] to place the midpoints of the frac- 



The outcrops can be divided into two groups: one with 



tures on a fractal set characterized by a fractal dimension 
D2- We show in Fig. [5] the resulting fracture systems for 
different parameters ai and -D2. The outcrop data has 
D2 ~ 2. Further details may be found in |22i] . 

The results of analyzing the equivalent networks of the 
DFN model networks are given in Table [Hi The data are 
based on 1000 networks of comparable size to those in the 
outcrop fracture data sets. From the table, we see the 
same trends as those observed in Table U for the eight 
outcrop fracture data sets and it is possible to find a 
combination of a/ and D2 to make match between them. 
However, we show in Fig. [5] the averaged degree cor- 
relation matrix. This indicates an assortative network 
structure: nodes of equal coordination number tend to 
be connected. This is the opposite of what is observed 
for the outcrop data sets, see Fig.H) Hence, the topology 
of the artifical networks is quite different from the natu- 
ral ones. This implies that the topology of the fracture 
network themselves, artifical and real, are quite different. 
This difference is not visible from direct observation. 

Hence, by constructing the equivalent networks, we 
have access to the entire analysis toolbox of modern net- 
work theory for fracture networks. As we have shown 
in the analysis presented, this makes it possible to test 
fracture network models on a quantitative level beyond 
what has been possible earlier. 

We thank Svensk Karnbranslchantering AB for out- 
crop data, H. F. Hansen and E. Skjctne for discussions 
and the referees for valuable comments. C. A. A. and A. 
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FIG. 6: The degree correlation matrix C{ki, ^2) for different 
DFN model parameters a; and D2. This figure should be 
compared with Fig. |4l 



TABLE II: List of degree distribution power-law exponents 
ak, clustering coefficient C, clustering coefficient for rewired 
networks Crw, clustering coefficient for comparable random 
networks Cra, efficiency E, efficiency for rewired networks 
Erw, and efficiency for comparable random networks Era 
for various fracture length power-law exponents ai. 

ai ctk C Crw Cra E Erw Era 

2.00 2.2 0.08 0.0f9 0.047 0.028 0.042 O.lf 
2.25 f.7 O.lf 0.013 0.031 0.027 0.049 0.11 
2.50 1.4 0.17 0.013 0.019 0.037 0.083 0.10 
2.75 1.2 0.26 0.014 0.014 0.050 0.134 0.09 
3.00 0.9/1.3° 0.31 0.013 0.008 0.050 0.154 0.07 

"A kink in the slope around fc = 60 gives 0.9 when fitting for 
smaller values of k and 1.3 for larger values. 
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